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ABSTRACT 

We propose a method to infer the initial spin period of pulsars residing in composite supernova remnants. 
Such a remnant consists of both a plerionic and a shell type component, corresponding respectively to the pulsar 
wind nebula driven by the spindown luminosity of the central pulsar, and the blastwave bounding the supernova 
remnant. Theoretical investigations including hydrodynamical simulations have shown that at late times 1,000 
- 10,000 years), a simple scaling law connects the radius of the supernova shell to the radius of the plerion. The 
energy content of the plerion and the total mechanical energy of the supernova remnant enter into this scaling law. 
One can use this scaling law to estimate the initial spin period of pulsars residing in composite remnants. We 
discuss potential pitfalls of this method, including the effect of a small remnant age and of strong radiative losses 
in the plerion. 

Subject headings: hydrodynamics — (ISM:) supernova remnants — (stars:) pulsars: general 
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1. INTRODUCTION 

Is a typical neutron star born spinning close to break-up 
(P() ~ 0.8 ms, Haensel et al. 1995). The answer to this ques- 
tion reveals much about the angular momentum evolution of 
neutron star progenitors (Heger et al. 2000), as well as about 
the process es tha t occur during the super nova explosion (Spruit 
& Phinney 119981 Lai & Goldreich|2000|). If the initial spin-rate 



is close to the break-up rate, the nascent neutron star may un- 
dergo a Rossby-wave instability (Lindblom, Owen & Morsink 
1998) and emit gravitational radiation. Moreover, a rapidly ro- 



tating neutron star can obtain a substantial kick velocity through 
the electromagnetic rocket effect (Harrison & Tademaru 1975). 



This can provide a natural explanation for the apparent align- 
m ent be tween the spin axis and pulsar pr oper m otion (Lai et 
al. [200 1[ ) in the Crab (C araveo & Mignani |1999[ ) and Vela pul- 
sar (Pavlov et al. [2000| ). 

A number of studies have been dedicated to the determina- 
tion of initial spin-rates. Those few young pulsars which have 
measured braking indices and estimated ages s eem t o possess 
initial spin periods from 19 ms (Crab, Lyne et al. 1993[ ) to 63 ms 
(PSR B 1509-58, Kaspi et al. |1994| ).^ Since only rapidly rotating 
pulsars have measurable braking indicies, these values may not 
be representative f or the general population. For these, Phin- 



ney & Blandford (1981) and Vivekanand & Narayan (1981) 



have developed the pulsar-current analysis. Unfortunately, 
thi s met hod is subject to small-number statistics (Lorimer et 
al. 1993). Analysis based on pulsar luminosity functions (Em- 
mering & Chevalier 1989, Narayan |l 987 ) raised the possibility 
of 'injection' : the majority of the pulsars may be born with peri- 
ods as slow as several hundred milliseconds. This hypothesis is 
corroborated by the ap paren t paucity of plerions around neutron 
stars (Srinivasan et al. [1984 ). A plerion derives energy from the 



rotational energy of the central pulsar. The overall conclusion 
from pulsar studies seems to favour a slow initial spin rate. 

In this article, we propose a new method for inferring the 
initial spin-rates of some neutron stars - the small population 
that reside within composite supernova remnants. A compos- 



ite supernova remnant (SNR) includes both a plerionic (filled- 
centered) component (Weiler & Panagia 1978 ) and a shell com- 



ponent. The first corresponds to the pulsar wind nebula (PWN), 
while the latter corresponds to the blastwave of the supernova 
remnant (SNR) propagating through the interstellar medium 
(ISM). They are observationally distinguished by a different 
spectral power-law index at radio frequencies. The dynamics 
of such a composite system has been extensively studied, e.g. 
Reynolds & ChevaUer (1984); Chev alier & Fransson (1992). 

Recently, van der Swaluw et al. (2001) presented hydrody- 
namical simulations and analytical arguments which explicitly 
relate the radius of the PWN to that of the SNR. These authors 
find that the two radii are roughly proportional some time af- 
ter the initial explosion, with the proportionality constant deter- 
mined by the ratio of the pulsar spin-down energy and the total 
mechanical energy of the supernova event. This relation forms 
the theoretical basis for the current paper. In this article, we 
follow van der Swaluw et al. (2001 ) in restricting ourselves to 
spherically symmetric systems. 

2. EVOLUTION OF A PWN IN A SNR 

As a pulsar spins down, its rotational energy is largely de- 
posited into the surrounding medium, driving a relativistic pul- 
sar wind into the tenuous bubble of stellar debris left behind by 
the supernova blastwave. We approximate the spin-down lumi- 
nosity of the pulsar wind by that of a rotating magnetic dipole: 



L{t) = La/ (l+t/T) 



(1) 



The total spin-down energy is E^^ = Lqt. We neglect the dy- 
namical influence of the pulsar wind on the SNR, assuming that 
Esd is well below the total mechanical energy of the supernova 
explosion, Eq 10^' erg. 

The evolution of a PWN inside a SNR can be divided into 
two stages (Reynolds & Chevalier 1984; Chevalier & Frans- 
son 1992). In the first stage, the PWN expands supersonically 
within the bubble blown out by the SNR. The second stage 
commences at f = tsj when the PWN encounters the reverse 
shock of the SNR. The reverse shock heats the interior of the 
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SNR and causes the expansion of the PWN to become sub- 
sonic. In this stage the expansion of the PWN is regulated by 
the expansion of the SNR, which is described by a self-similar 
Sedov-Taylor solution. In the following we summarize some 
relevant result s from Chevalier & Fransson (1992) and van der 
Swaluw et al. fcOOlh. 



2.1. Analytical Relations 

We focus on the subsonic expansion stage of the PWN which 
occurs when the SNR has relaxed to the Sedov solution. We 
assume that the pulsar wind has deposited most of its energy 
into the PWN (f > r). We find that the radius of the PWN 
(Rpv/n) scales roughly linearly with the radius of the SNR (^sm )- 
This scaling arises from the condition of pressure equillibrium 
between the PWN interior and the interior of the SNR. We ap- 
proximate the pressure which confines the PWN by t he cen tral 
pressure of the SNR, using the Sedov solution (Sedov 1958), 
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snr ' 



(2) 



where Eq is the total mechanical energy of the SNR. Inside the 
PWN, pressure is quickly equilibrated due to the high sound 
speed ~ c/-\/3 of the relativistic fluid. The interior pressure is 
related to the energy content of the PWN, by 



~ 1/3 

unil\ . The second parameter is defined as 773(f) = Cryj < C. 
The maximum value of 773 = C is used when we determine the 
initial spin periods of pulsars driving composite remnants. 

The derivation leading up to equation (jj) is invalid during 
the earlier supersonic evolution. We retain equation (Bh with 
the sole purpose of defining 773 during these stages. We pro- 
ceed to demonstrate the behaviour of 773 throughout the PWN 
evolution with results from hydrodynamical simulations. 



2.2. Hydrodynamical Simulations 



Van der Swaluw et al. ( |2001| ) have presented hydrodynam- 
ical simulations of a pulsar wind nebula inside a supernova 
remnant. In these simulations, the gasdynamical equations are 
integrated in a spherically symmetric configuration, using the 
Versatile Advection Code (VAC) developed by Gabor Tothat 
the Astronomical Institute Utrecht (Toth and Odstrcil 1996 ). 
Lacking the possibility of treating a relativistic fluid, the pul- 
sar wind has been implemented as a cold non-relativistic wind 
(7 = 5/3) with a terminal velocity equal to the speed of light, 
Voo = [2L(t)/Mpv/(t)V^^ — c, where Mp„(f) is the mass ejection 
rate into the pulsar wind bubble, and L{t) is the pulsar spin- 
down luminosity as given by equation (|l]). 
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Let Epwn be the total amount 



where 7 is the adiabatic index 
of energy injected into the PWN by the central pulsar. Part of 
this energy is used to perform work on the surrounding medium 
as the PWN expands. Despite some compression of the PWN 
by the reverse shock around t = f sr, -E* falls below Ep^/n (see 
Fig. 9 of van der Swaluw et al. ( 2001 )). We therefore write 
£■* = 772£pwn with 772 < 1 decreasing over time. Imposing pres- 
sure equilibrium, we find 



^pwn = C (E^/Eq) 



1/3 



(4) 



1.29 



where C ~ 1.02 for a relativistic fluid (7 = 4/3) and C : 
for a non-relativistic fluid (7 = 5/3). 

In the absence of radiative losses, the subsonic expansion of 
the PWN is adiabatic, PpwnVpwn = constant, where Vp^a « Rp^„ 
is the volume of the PWN. Combining this relation with equa- 
tion we find that falls off with time as for a rela- 
tivistic fluid. This yields an expansion law for the radius of the 
PWN, Rp^^ cx f^/"', wher eas the radius of the SNR scales as in 

a Sedov solution (/?sm t'^l^). So, roughly, Rp^^n oc /?sm- 

We introduce two additional dimensionless parameters, 771 
and 773, and rewrite equation (^ into a form that can be more 
easily compared with simulations. 



(5) 



The parameter 771 is a constant and relates the total energy input 
into the PWN, Ep„n, to the total spin-down energy of the pulsar, 
iisd as Ep„n = ViEsA- This takes into account possible radiative 
losses in the PWN, and other inefficiencies in the conversion of 
the spin-down energy into mechanical energy. For instance, a 
fraction of the neutron star spin-down energy may escape di- 
rectly from the pulsar as high energy radiation. We set 771 = 1 
in this article, but discuss the case when it is much less than 




^ 3 o 



Time (in kyears) 

Fig. 1 . — Results of numerical simulation showing the time evolution of the 
pulsar wind nebula. Plotted here are its radius (in units of parsec, thin contin- 
ued line) and the dimensionless parameter TI3 (thick line). The latter is disected 
into three parts: a supersonic expansion stage when the PWN is bounded by a 
shock (dotted section), in unsteady transition when undergoing reverberation 
with the supernova reverse shock (dashed section), and subsonic expansion 
(solid section). The parameters adopted in this simulation is: maximum spin- 
down luminosity Lq = 5 X lO'** ergs/s, spin-down time r = 600 yrs, supernova 
explosion energy of £"0 = 10^' srgs 2nd ISM density of 10"^* g/cm^. We find 
?ST ^ lO'yrs in this case. Note that rj} is not physical in the early expansion 
stages but is opl^iJprmally defined by equation (ph. For details, see van der 
Swaluw et al. ( pOOl| ). 

Figure |l] presents relevant results from one of these simu- 
lations. The numerical parameters adopted in this simulation 
are listed in the caption. Taking different parameters will not 
qualitatively change the overall behaviour of the system, but 
it will affect the early (supersonic) evolution and the moment 
when subsonic expansion commences (fsx)- For example, fsx 
is larger if the supernova explosion is more energetic, if the su- 
pernova ejecta mass is larger, or if the interstellar medium den- 
sity is smaller (fsj oc Eq^^^M^I^ pjj^'^ , e.g. McKee & Truelove 
(1995)). Since this is a non-relativistic simulation (7 = 5 /3), we 
expect ri2 <C = 1.29 in the subsonic expansion stage. This is 
indeed observed. Therefore in the following investigation we 
assume a maximum value for 773 of 1.02, as is appropriate for 
the subsonic expansion stage of a PWN in a SNR with 7 = 4/3. 



Inferring neutron star initial spin periods 
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3. INFERRING INITIAL SPIN RATES 

Writing the spin-down energy as Esa = (fll-ilj)I /2, we ob- 
tain the following expression for the pulsar's initial spin period. 



(6) 




where the spin-periods are Pq = 2tt/Q,(), = 2ti/Q,, and / is the 
moment of inertia of the neutron star The radius of a PWN rel- 
ative to its associated SNR shell can be used to infer the initial 
spin-period of the pulsar, assuming ?7i, 773, £0 and il, are known 
or can be estimated. 



1000 




Fig 2. — Initial spin periods inferred from the ratio Rp„„/R siir tising equa- 
tion m). When assuming no radiative loss (771 = 1), we obtain solutions as 
depicted by the family of solid lines, with different lines applicable for differ- 
ent cuiTent spin periods, Pi. The family of dashed lines are for 771 = 0.05. Filled 
and open triangles represent systems listed in Table M for »?i = 1 and rji = 0.05, 
respectively. When Pq <A ms, the initial rotationaT energy oLlhe puslar ex- 
ceeds the total mechanical energy in the SNR itself. Equation dq) fails and the 
SNR may be blown away by the PWN. None of the systems we examined lies 
close to this limit. 

The radius of PWN should be taken from radio measure- 
ments. Synchrotron-radiating electrons in PWN have a typical 
Ufe-time (Rybicki & Lightman[l979|): 



^1/2 = 



WB 



,-3/2 -1/2 

mC ^MHz yi"^' 



(7) 



£■() = 1.0 X 10^' erg, / = 1.4 x lO'*^ g cm^ (the value for Crab), 
and P, = 211/0,, = 00 in cases where it is unknown. The same 
results appear in Figure || as filled triangles. By taking the max- 
imum value for 773 we probably under-estimate the initial spin 
rate, but likely by no more than a factor of 2 judging from Fig- 
ure |l} In the same table we also list the initial spin periods {Pq) 
derived for 4 pulsars using braking index measurements. 

4. DISCUSSION 

We have demonstrated a new method for inferring initial spin 
rates of pulsars residing in composite supernova remnants. This 
method uses the ratio of the plerion radius and the radius of the 
supernova shell, so we dub it as the 'radius method' . Given that 
typical plerions and supernova remnants are expected to live for 
lO-'- 10"* years, we expect that future observations will increase 
the population to which the radius method can be applied. 

This method does not require a knowledge of parameters like 
the density and magnetic field strength of the ISM, or the dis- 
tance and age of the SNR. However, for the method to apply, 
the observed PWN has to be in the subsonic stage of its evo- 
lution. If this is not the case, the method systematically over- 
estimates the initial spin rate (see below). Uncertainties in Eq, 
rj-i and VL, introduce errors of order unity, as will a deviation 
from pure spherical symmetry. A more important uncertainty 
of this new method concerns the influence of synchrotron ra- 
diative losses in the plerion, as parametrized by r/i . We tend 
to under-estimate the initial spin-rate by ignoring these losses 
(setting rii = 1). We will illustrate the effect of these uncertain- 
ties by considering the 4 systems in Table [l|, which have known 
central pulsars, and independent estimates for their initial peri- 
ods based on measurements of the pulsar breaking index. ^ 

Among these systems. Vela is the oldest and should have en- 
tered the subsonic expansion stage. Its associated pulsar being 
displaced from the center of the plerion also testifies in favor of 
an old age. Adopting r/3 = 1 .02 and other parameters as stated, 
the radius method yields an initial spin-down energy fo r Vela 
that is a factor of 2 above that inferred by Lyne et al. (1996) 



where ctj is the Thomson cross-section, F is the electron 
Lorentz factor, v is the typical synchrotron frequency scaled 
at the radio wavelength, and B is the magnetic field strength 
scaled by a value appropriate for plerions. Such a long life-time 
implies that radio observations 'see' electrons that were accel- 
erated during the early energetic stage when the pulsar deposits 
most of its rotational energy into the pulsar wind. They there- 
fore delineates the true spatial extent of the plerion. In contrast. 
X-ray electrons loose their energy quickly. The X-ray plerion 
is produced by electrons recently accelerated near the pulsar, so 
one ecpects the plerion to be smaller in X-rays. 

We have collected from the literature 13 composite super- 
nova remnants. We list the relevant properties of these sys- 
tems in Table |l|. The relative sizes of the plerions are taken 
from radio observations, with the exception of Gil .2 - 0.3 and 
G320.4- 1.2 (marked by asterisks), where only X-ray plerions 
have been observed. We apply equation (||) to determine the ini- 
tial spin period for these systems, adopting ryi = 1 .0, 773 = 1 .02, 

^GI 1.2— 0.3 does not have a braking index measurement. But its inferred initial 



from pulsar braking index measurement. This is surprising as 
we have ignored radiative loss. We suggest three possibilities 
that will reconcile the two methods: (1) the supernova explo- 
sion is half as energetic as assumed here, (2) the pulsar is 15 
ky r old, and (3) the braking index is as high as 2.4 (Lyne et 
al. 1996 reported 1.4). The first option will raise Pq to 52 ms 



while the latter two will reduce P'^ to 41 ms. 

The remaining two systems, Gil. 2 -0.3 and G320.4- 1.2, 
are young (^ 1.5 kyr). Their plerions could still be in super- 
sonic expansion during which the value of 773 (defined formally 
through equation (^) in the supersonic stage) can be higher than 
1 .02, the theoretical maximum during the subsonic expansion 
stage (Fig. |l]and §2 of van der Swaluw et al. 2001). For exam- 
ple, taking 773 = 2 we find Pq = 63 ms for Gl 1.2-0.3. A similar 
correction may have to be applied to G29.7-0.3 (Kes 75). 

In contrast to the systems discussed above, SNR 0540-69.3 
is the only case in which the spin-down energy is significantly 
under-estimated for our choice of parameters. The plerion in 
SNR 0540-69.3 is too small if the pulsar w as inde ed born with 
an initial spin of 39 ms. Manchester et al. (1993) have argued 



that the radiative loss in this system has reduced the total energy 
stored in the plerionic magnetic field and relativistic electrons 
to a mere 4% of the integrated spin-down luminosity of the pul- 
sar. Taking rji = 0.04 we recover a 39 ms initial period, with a 

spin depends little on this index due to its small age. 
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Supernova 
Designation 



Plerion 



Pulsar 



Pa (in msec) 



Pt (in msec) 



Pq (in msec) 



age 
(kyr) 



GO.9+0.1 

G16.7 + 0.1 

G29.7-0.3 

G34.7-0.4 

G293.8+0.6 

G322.5-0.1 

G326.3-1.8 

G327.1-1.1 

G351.2+0.1 



0.25(Helfand & B ecker 



1987) 



0.25(Helfand et al. |1989E 
0.16(Becker& Helfa nd 
0.07 (Frail etal.|l99" 



1984) 



0.25(Whiteoak & Gr een |1996| ) 
0.27(Whiteoak |l99ji 
0.28(Kassim et aT 
0.25(Sun et al. 



1993) 



19^) 

0.05(Becker & Helfand|l988|) 



43 
43 
82 

197 
43 
39 
37 
43 

484 



325[PSR 1846-0258] 
267[PSR 1853+01] 




0.7 



G263.9-3.3 
Gll.2-0.3* 
G320.4-1.2* 
SNR0540-69.3 



0.24(Milne 1980) 



0.15(Vasisht et al. 1996 i 
0.17(Seward et al. H983|) 
0.08(Manchester etal.|1993|) 



41 
53 
69 
49 



89 [Vela] 
65 

150[PSRB 1509-58] 
50.3[PSR 0540-69] 



52(Lyne et al. 
63(Torii et al. 
63(Kaspi et al. 
39(Deeter et al 



11.3 
1.6 
1.5 
0.8 



Table 1 

List of composite remnants: Po is the initial spin period from the radius method, Pt is the current spin period of 

SYSTEMS WITH IDENTIFIED PULSARSQ AND Pq IS THE CALCULATED INITIAL SPIN PERIOD BY USING THE BRAKING INDEX. 



value of 773 = 1.3. This value of 773 is reasonable because SNR 
0540-69.3 is only 760 years old. None of the other 3 sys- 
tems discussed here requires strong synchrotron losses in the 
plerion, which makes SNR 0540-69.3 unique. Among the 4 
pulsars, PSR 0540-69 is not outstanding either in its surface 
magnetic field strength or in its spin-down rate. It is likely born 
the fastest, though not by a larg e marg in. 

Along the same line, Atoy an ( 1 999 ) have argued that in Crab, 
strong synchrotron losses in its past can account for the ob- 
served flat spectra of the radio electrons in the plerion. If this 
appUes to all PWNe, it would significantly change our infeiTed 
initial spin rates. This is illustrated in Fig. |^ for a strong syn- 
chrotron loss of ri\ = 0.05. 

Although the effect of radiative losses on the reliability of our 
results remain unclear, we briefly explore the implications of 



ny 

these results. Our sample (table nj) is likely biased towards fast 



spinners, which produce plerions that are easier to detect. Even 
so, all pulsars in our sample seem to be born spinning much 
below the break-up rate. Moreover, their initial spin-periods do 
not cluster around 5- 10 ms, a signature of rapid spin-down 
by the gravitationally radiating Rossby-waves which occ ur im- 
mediatel y after the supernova collapse (Andersson et al. 2000, 
Wu et al. 2001). The long initial periods we find imply an ef- 
fective angular momentum coupling between the core and en- 
velope of neutron star progenitors until late into their evolution. 
The scattering in the initial periods could either reflect a differ- 
ent decoupling time in different stars, or a stochastic process 
that gives rise to the final angular momentum in the neutron 
star, suc h as an off-centered kick during the collapse (Spruit & 
Phinney[T998l). 



We acknowledge discussions with Bram Achterberg, Armen 
Atoyan and an anonymous referee. 
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